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Abstract 

Solvent  effects  on  the  rate  constants  for  both  homogeneous  and  heterogeneous 
electron  transfer  reactions  have  been  analyzed  on  the  basis  of  current  models  which 
consider  the  role  of  dynamic  relaxation  processes  in  determining  the  magnitude  of  the  pre¬ 
exponential  factor.  A  statistical  method  for  separating  the  effects  of  the  solvent 
longitudinal  relaxation  time  tl  from  those  of  the  solvent  permittivity  parameter  is  ‘ 

described  and  applied  to  15  sets  of  experimental  data  for  which  results  are  available  in  at 
least  four  solvents.  The  degree  to  which  the  explained  variation  in  the  logarithm  of  the  rate 
constant  could  be  attributed  to  either  of  these  effects  varied  all  the  way  from  0  to  100% 
depending  on  the  degree  of  reaction  adiabaticity  and  the  relative  sizes  of  the  inner  and  outer 
sphere  components  of  the  free  energy  of  activation.  Data  for  the  limiting  cases  in  which 
there  is  no  Xl  dependence  in  the  pre-exponential  factor  or  in  which  the  pre-exponential 
factor  is  proportional  to  were  analyzed  further  to  obtain  the  size-distance  parameter  and 
the  components  of  the  pre-exponential  factor  relevant  to  the  encounter  pre-equilibrium 
model.  These  parameters  have  been  discussed  with  respect  to  current  developments  in 
electron  transfer  theory.  Problems  in  estimating  the  longitudinal  relaxation  time  in  the 
solvent,  required  for  the  analysis,  are  also  considered. 

y\ 
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Introduction 

The  importance  of  solvent  dynamical  properties  in  electron  transfer  reactions  has  been 
recognized  and  elaborated  in  recent  theoretical  work  [1-5].  At  the  same  time,  solvent 
effects  on  the  kinetic  parameters  for  simple  electron  transfer  reactions  involving  molecular 
reactants  and  their  corresponding  monovalent  cations  or  anions  have  been  studied  for  both 
homogeneous  [6- 10]  and  heterogeneous  processes  [1 1-19]  in  a  wide  variety  of  solvents. 
Such  systems  are  ideal  for  investigating  solvent  effects  because  work  terms  are  a 
minimum,  and  the  inner  sphere  free  energy  of  reorganization  is  often  much  smaller  than  the 
outer  sphere  quantity.  In  earlier  work,  it  was  assumed  that  the  solvent's  primary  role  was 
an  activational  one,  as  expressed  in  the  now  famous  Marcus  dielectric  continuum  treatment 
[20].  Certainly,  the  first  studies  that  considered  a  variety  of  solvents  for  homogeneous 
[6]  and  heterogeneous  [1 1,21]  processes  showed  that  the  early  form  of  the  equations  for 
the  rate  constant  in  which  the  pre-exponential  factor  is  not  solvent  dependent  failed 
dramatically.  This  was  in  stark  contrast  with  the  results  of  intramolecular  electron  transfer 
processes  for  which  excellent  agreement  was  shown  in  many  cases  between  charge-transfer 
band  adsorption  maxima,  and  the  quantity  predicted  from  electron  transfer  theory  [22]. 
These  studies  offer  very  strong  evidence  for  the  general  validity  of  the  Marcus  treatment 
[20]  in  so  far  as  the  activation  process  is  concerned.  Inspired  by  the  recent  theoretical 
work,  more  effort  has  been  expended  in  understanding  the  role  of  the  solvent  in 
determining  the  pre-exponential  factor  [7-10,12,14-19]  in  the  rate  constant  for 
intermolecular  electron  transfer.  The  evidence  is  now  quite  clear  that  both  homogeneous 
and  heterogeneous  kinetic  data  correlate  quite  well  with  solvent  relaxation  times 
[15,17,23].  However,  the  precise  functional  dependence  of  the  electron  transfer  rate 
constant  on  solvent  relaxation  time  is  by  no  means  clear  at  the  present  time.  Early  work  by 
Zusman  [1]  predicted  that  the  pre-exponential  factor  in  the  rate  expression  is  proportional  to 
Xl**  where  Xl  is  the  longitudinal  relaxation  time.  Recent  work  by  Marcus  [5,24]  has 
shown  that  Zusman' s  formulation  is  a  limiting  case,  such  that,  as  the  inner  sphere 
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reorganizadonal  energy  increases,  the  nuclear  frequency  factor  Vn  should  acquire  a 
fractional  dependence  on  that  is,  Vn  is  proportional  to  Xl‘“  where  a  is  a  fraction 
between  0  and  1.  In  addition,  if  the  reaction  is  weakly  adiabatic,  one  also  expects  the 
dependence  of  the  pre-exponential  factor  on  Xl'*  to  be  weakened  such  that  for  non-adiabatic 
processes  the  solvent  dependence  of  the  pre-exponential  factor  disappears  [25,26]. 

With  much  of  the  attention  focussed  on  solvent  dynamical  effects,  we  were  interested 
in  determining  whether  the  Marcus  dielectric  continuum  model  for  the  outer  sphere 
reorganizadonal  energy  can  be  tested  for  simple  electron  transfer  reactions  using  kinetic 
data  obtained  in  a  number  of  solvents.  In  a  recent  paper  [27],  we  analyzed  data  for  solvent 
effects  on  electrochemical  rate  constants  on  the  basis  of  the  limiting  form  of  the  solvent 
dependence  of  the  pre-exponential  factor  (a  =  1),  and  found  that  the  reorganizadonal  free 
energy  of  activation  depends  on  the  solvent’s  dielectric  properties  in  a  manner  predicted  by 
electron  transfer  theory  [20].  As  a  result  of  our  analysis,  it  was  possible  to  assess  other 

parameters  relevant  to  electron  transfer  theory  such  as  the  size-distance  parameter  which 
determines  the  magnitude  of  the  activation  energy,  and  the  parameters  for  precursor 
complex  formation.  Another  important  test  of  electron  transfer  theory  involves  a 
comparison  of  kinetic  data  for  homogeneous  and  heterogeneous  electron  transfer.  Such  a 
comparison  necessarily  involves  assessing  the  role  of  imaging  in  determining  the 
magnitude  of  the  size-distance  parameter  for  the  heterogeneous  process.  One  also  needs  a 
way  of  assessing  the  importance  of  solvent  dynamical  effects  in  determining  the  pre¬ 
exponential  factor.  In  the  present  paper,  we  extend  our  previous  analysis  [27]  to 
determine  the  importance  of  solvent  dynamical  effects  in  the  kinetic  data  and  compare 
results  for  homogeneous  and  heterogeneous  reactions  using  data  presented  in  the  literature. 
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Theory  and  Analysis 

According  to  Marcus'  theory  for  electron  transfer  reactions  [20,28,29],  the  rate 
constant  for  electron  transfer,  kf,  may  be  written 

kr  =  kZ  exp  (-AG*/RT)  (1) 

where  ic  is  the  electronic  transmission  coefficient,  Z,  a  collision  frequency  and  AG*,  the 
reorganizational  free  energy  of  activation.  It  is  understood,  in  the  present  case,  that  the 
term  Z  contains  any  electrostatic  work  terms  associated  with  bringing  the  reactant  or 

reactants  to  the  reactant  site,  so  that  the  activation  free  energy  is  composed  of  two  terms, 

*  * 

AGjj,  the  inner  sphere  contribution  due  to  internal  structural  reorganization,  and  AG^^,  the 

* 

outer  sphere  contribution  from  the  solvent  medium.  AG-^  can  be  estimated  from  the  bond 

coordinate  changes  between  the  oxidized  and  reduced  forms  of  the  reactant,  and  the 

corresponding  bond  force  constants  which  can  be  determined  firom  vibrational  spectroscopy 

* 

[28,29] .  The  outer  sphere  quantity  AG^^  for  homogeneous  bimolecular  reactions  is  given 
by  [20,29] 


where  No  is  Avogadro's  number,  e,  the  electronic  charge,  Eq,  the  vacuum  permittivity, 

ai,  and  a2,  the  radii  of  the  two  reactants  represented  as  spheres,  Rh,  their  separation  in  the 

encounter  complex,  and  Eop  and  £5,  the  optical  and  static  dielectric  constants  of  the 

solvent,  respectively.  Since  we  are  concerned  here  with  reactions  involving  a  molecule, 

and  its  corresponding  cation  or  anion,  ai  =  82  =  a,  and  the  size-distance  parameter  in  eq. 

* 

(2)  reduces  to  (a*^  -  Rh'^).  The  corresponding  value  of  AG^^  for  heterogeneous  or 
electrochemical  reactions  is 


327t£o  leop  ’  EsJ 


(3) 


where  Re  is  the  distance  between  the  reactant  and  its  image  in  the  conducting  electrode. 


6 


These  expressions  for  are  based  on  the  Bom  model  for  ion  solvation,  and  hence, 

inherit  some  key  assumptions,  namely,  that  the  solvent  is  a  structureless  continuum  and 
that  there  is  no  spatial  dispersion  of  the  dielectric  permittivity  at  distances  beyond  the 
reactant  radius. 

In  its  original  form,  the  Marcus  expression  [20]  for  the  homogeneous  rate  constant 
treated  the  pre-exponential  factor  in  a  'gas  phase'  manner  as  an  effective  collision 
frequency,  Zj,.  When  one  includes  the  work  term,  the  expression  for  Zh  may  be  written 

D'T  1/2 

Zh  =  4:cNoRhY - exp(-wh/RT)  (4) 

\^27cM  j 

where  wh  is  the  electrostatic  work  done  to  bring  the  reactants  to  the  reaction  site 

and  M  is  the  reduced  mass  of  the  reactants.  In  the  case  of  heterogeneous  processes,  the 

collision  frequency  was  written  as 

DT  1/2 

^  exp(-We/RT)  (5) 

where  We  is  the  electrostatic  work  done  to  bring  the  reactant  to  the  reaction  site  near  the 
electrode  and  M  is  now  understood  to  be  the  mass  of  the  reactant. 

On  the  basis  of  the  more  recent  encounter-preequilibrium  model  [25],  the  pre¬ 
exponential  factor  in  equation  (1)  has  been  recast  as  KKpVp  where  Kp  is  the  equilibrium 
constant  for  precursor  complex  formation  and  Vn,  the  nuclear  frequency  factor  giving  the 
effective  frequency  with  which  the  system  crosses  the  free  energy  barrier.  In  the  case  of 
homogeneous  electron  transfer  [25],  the  equilibrium  constant  Kp  is  given  by 

Kp  =  47C  NoRh^  5r  exp  (-Wh/RT)  (6) 

where  5r  is  the  range  of  reactant  separations  over  which  electronic  coupling  is  sufficient  for 


reaction.  The  corresponding  equation  for  a  heterogeneous  process  [30]  is 

Kp  =  5r  exp  (-We/RT) 


(7) 
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The  magnitude  of  5r  which  is  important  in  determining  Kp  for  both  cases  is  considered  to 
vary  between  0.2  nm  for  adiabatic  reactions  to  30  pm  for  non-adiabatic  reactions  [31].  It 
should  be  noted  that  the  formulation  for  the  pre-exponential  factor  used  here  does  not 
explicidy  include  the  nuclear  tunnelling  factor  [25,31]  which  can  be  assumed  equal  to  unity 
for  most  electron  transfer  reactions  at  room  temperature. 

In  general,  the  nuclear  frequency  factor  Vn  depends  on  both  internal  vibrational  and 
external  solvent  flucmation  frequencies  for  the  reactant(s).  A  convenient  way  of 
expressing  Vn  is  in  terms  of  a  weighted  average  of  frequencies  associated  with  inner  sphere 
vibrations,  Vis  and  outer  sphere  solvent  reorientation,  Vqs,  that  is  [25,31], 


where  eoo  is  the  'infinite  frequency'  dielectric  constant.  The  quantity  £«  is  regarded 
differendy  by  different  authors,  some  [5,33]  defining  it  as  the  optical  dielectric  constant 
Eop  (n^  where  n  is  the  refractive  index)  and  others  [16-18]  as  the  value  in  the  infrared 
frequency  range.  This  problem  will  be  addressed  further  below.  Finally,  in  the  limit  that 

4i  >)■ 

AGjg  is  much  greater  than  AG^^,  Vn  is  equal  to  Vis  and  the  solvent  dependence  of  the  pre¬ 
exponential  factor  disappears  [5,24]. 
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A  criterion  for  determining  whether  solvent  dynamical  frequencies  should  predominate 
in  the  pre-exponential  factor  of  eq.  (1)  was  given  by  Ovchinnikova  [3,13] ;  accordingly  , 
eq.  (9)  is  valid  if  the  following  condition  is  met: 


Vis  [AG*/ AG*)  +  AG*s)]1/2  exp(-AG*/RT)  >  Tl’I 


(11) 


Certainly,  many  systems  are  not  expected  to  fulfill  this  criterion.  Recent  work  by  Nadler 
and  Marcus  [24]  predicts  that,  as  AG.^  increases  relative  to  AG^^,  the  rate  constant  kj. 
assumes  a  power  law  dependence  on  Xl  such  that  it  is  proportional  to  Xl'“  where  a  is  a 
fraction  between  0  and  1.  On  the  basis  of  this  work,  a  general  way  of  writing  eq.  (1)  for 
adiabatic  reactions  with  inclusion  of  solvent  dynamical  effects  in  the  pre-exponential  factor 
is 


kr  =  A  Xl-“  exp  (-AG*/RT)  (12) 

*  * 

where  the  parameter  A  is  equal  to  icKp  Vjs  when  AGj^  »  AG^^,  and  to  xKpVos  when 
*  * 

AGqj  »  AGjj,  and  the  parameter  a  is  a  fraction  between  0  and  1.  The  same  expression 
may  be  used  for  weakly  adiabatic  reactions  for  which  the  coefficient  a  is  expected  to 


decrease  from  unity  with  departure  from  adiabaticity  [26].  Rewriting  eq.  (12)  making  use 


of  the  expressions  for  AG^^  given  above,  one  obtains  the  result 


AG”^ 


In  kr  =  In  A  -  -  a  In  xl  -  g  Y 


(13) 


In  this  equation,  y  is  the  permittivity  parameter  ’tiven  by 


and  the  quantity  'g'  is  defined  for  homogeneous  reactions  by 


(14) 


Noe2  _  j_-v 

~  ItittEoRT  /  ’ 


or  for  heterogeneous  reactions  by 

Noe2 

^  327teoRT 


(15) 


(16) 
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It  follows  that  one  may  assess  the  relative  importance  of  solvent  effects  in  the  pre- 
exponential  factor  and  in  the  outer  sphere  free  energy  barrier  by  performing  a  two 
parameter  least  squares  analysis  of  the  dependence  of  Inkf  on  IntL  and  y.  Strictly 
speaking,  the  quantity  InA  depends  on  y  for  adiabatic  reactions  with  a  =  1;  however,  this 
dependence  correlates  well  with  y  as  will  be  shown  below,  so  that  the  analysis  of  solvent 
effects  on  the  basis  of  eq.  (13)  provides  an  excellent  method  of  assessing  the  relative 
imponance  of  the  two  contributions.  It  should  also  be  kept  in  mind  that  the  parameter  a 
can  change  with  solvent  nature  for  a  given  reaction  [24].  Thus,  when  a  is  determined  by 
a  least  squares  fit  of  eq.  (13)  to  experimental  data,  the  result  obtained  should  be  regarded 
as  an  average  value  especially  when  it  is  not  at  either  of  the  limiting  values  (0  or  1). 

In  the  limit  that  inner  sphere  reorganization  energy  dominates  (a  =  0),  the  expression 


for  the  rate  constant  becomes 

In  kf  =  In  K  KpVis 


gY 


(17) 


Thus,  the  logarithm  of  the  rate  constant  should  be  linear  in  the  permittivity  parameter  y,  the 

size  distance  parameter  being  obtained  from  the  slope,  and  the  pre-exponential  factor  from 

* 

the  intercept  provided  the  value  of  is  known.  Such  a  plot  was  used  by  Grampp  and 
Jaenicke  [7]  to  analyze  the  solvent  dependence  of  kinetic  data  for  electron  transfer  involving 
diaminobenzenes  and  their  cation  radicals.  In  the  case  that  the  outer  sphere  reorganization 
energy  dominates  and  a  =  1,  it  can  be  shown  [23,27]  using  eqs.  (9)  and  (12)  that 


It  follows  that  a  plot  of  In  (kfXL/y^^)  against  y  should  be  linear  with  a  slope  of  -g. 

♦ 

Furthermore,  if  AGj^  can  be  estimated  from  spectroscopic  data,  one  may  estimate  icKp 
from  the  intercept  of  such  a  plot  [27].  It  should  be  pointed  out  that  eq.  (14)  assumes  that  if 
the  reaction  is  adiabatic,  the  degree  of  donor-acceptor  electronic  coupling  has  not  reduced 
the  barrier  height  or  altered  the  barrier  geometry  significantly  [4, 47]. 
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The  analysis  used  here  has  only  been  applied  to  data  obtained  in  aprotic  solvents. 

Most  protic  solvents  exhibit  multiple  relaxations,  and  there  is  some  theoretical  justification 
to  define  the  longitudinal  relaxation  time  in  terms  of  the  optical  dielectric  constant  rather 
than  the  infrared  value  [5].  Indeed,  Opallo  [17]  has  shown  clearly  a  difference  in  the  way 
heterogeneous  kinetic  data  for  aprotic  and  protic  solvents  correlate  with  the  value  of 

Obviously,  the  model  for  these  protic  solvents  must  be  more  complex  than  that  presented 
above,  and  therefore,  they  are  not  considered  in  the  present  paper.  Values  of  Xi  and  y  for 
the  solvents  involved  in  tne  following  analyses  are  given  in  Table  I.  The  sources  of  the 
original  data  and  values  of  Es,  eop>  £«.  and  Td  are  tabulated  elsewhere  [34].  For  the 
solvents  considered,  the  longitudinal  relaxation  time  varies  from  0.2  to  8.8  ps,  that  is,  by  a 
factor  of  50.  On  the  other  hand,  the  permittivity  parameter  varies  by  much  less,  and  falls 
within  the  range  0.28  to  0.53.  It  should  be  noted  that  homogeneous  electron  transfer  data 
[7,8]  have  been  reported  in  solvents  with  quite  small  values  of  the  permittivity  parameter  y. 
When  solvents  with  low  dielectric  constants  contain  an  electrolyte,  they  have  dielectric 
relaxation  characteristics  with  multiple  relaxations  due  to  ion  pair  formation  [35]. 

However,  the  data  reported  in  such  solvents  and  used  below  were  obtained  in  the  absence 
of  supponing  electrolyte  [7,8].  Finally,  the  analyses  represented  by  eq.  (13),  (17),  and 
(18)  should  strictly  speaking  be  applied  to  kinetic  data  which  have  been  corrected  for  the 
work  terms  introduced  here  in  the  equilibrium  constant  Kp  (eqs.  (6)  and  (7)).  The  work 
terms  are  normally  estimated  on  the  basis  of  the  extended  Debye-Huckel  theory  for  ionic 
strength  effects  in  the  case  of  homogeneous  processes  [25],  or  the  Gouy-Chapman  model 
for  double  layer  effects  in  the  case  of  heterogeneous  processes  [34].  Since  discussion 
below  for  homogeneous  reactions  is  limited  to  systems  in  which  one  of  the  reactants  is 
uncharged,  wj,  is  zero.  The  work  term  for  heterogeneous  electron  transfer,  Wg,  is  always 
non-zero  within  the  context  of  the  Frumkin  model  for  double  layer  effects  [36,37].  It  has 
been  argued  in  several  studies  [15-18]  that  double  layer  effects  are  negligible  for  the 
heterogeneous  reactions  considered.  This  subject  is  discussed  in  more  detail  below. 
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Results  and  Discussion 

The  dependence  of  the  rate  constant  for  electron  transfer  on  solvent  was  analyzed  on  the 
basis  of  data  for  both  homogeneous  and  heterogeneous  systems  which  have  already  appeared 
in  the  literature  using  the  equation 


Inkr  =  Inkro  -  alniL  -  gy  (19) 

In  this  equation  which  follows  directly  from  eq.  (13)  above,  the  parameter  kro  may  be 
regarded  as  the  solvent  independent  part  of  the  rate  constant  The  coefficients  a  ad  g  as  well 
as  the  normalized  partial  regression  coefficients  a  and  g  [38, 39]  were  determined  by  a  least 
squares  fit  of  the  experimental  data  using  the  values  of  Xl  and  y  summarized  in  Table  1. 

The  partial  regression  coefficients  a'  and  g'  are  defined  by  the  equations  [38,  39] 


and 


a'  =  lal 


<  In  XL  >) 


1/2 


<  In  kf  >) 


(20a) 


g'  =  Igl 


2(lnyi  -  <  y>)  \ 

j _ _ 


1/2 


\  ^  (In  kri  -  <lnk  >)/ 


(20b) 


1 

where  <ln  Xl>,  <y>,  and  <ln  kr>  are  average  values  of  the  quantities  In  Xl,  y,  and  In  kj, 
respectively,  and  the  summation  is  performed  over  the  parameters  involved  for  the  solvents 
used  in  a  given  set  of  experimental  data.  Since  a'  and  g'  are  on  the  same  scale,  comparison 
of  these  parameters  allows  one  to  assess  the  relative  importance  of  In  xl  and  y  in  determining 
the  observed  solvent  effect.  The  normalized  partial  regression  coefficients  defined  by  the 
equations 


a  = 


g' 

a'+  g’ 


(21a) 


and 


a'+  g' 


g  = 


(21b) 
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give  a  direct  estimate  of  the  fraction  of  the  variation  in  In  kr  with  solvent  explained  by  eq. 

(19)  which  can  be  attributed  to  variation  in  In  Tl  and  y,  respectively.  Although  the  actual 
value  of  a  is  important  in  assessing  a  given  set  of  data  and  is  reported  in  the  following  tables, 
very  little  significance  can  be  attached  to  the  value  of  g  obtained  in  this  analysis  especially 
when  a  is  close  to  or  equal  to  unity.  This  follows  from  the  fact  that  the  parameter  A  in  eq. 
(13)  depends  on  when  a  =  1,  and  that  there  is  a  strong  correlation  between  In  yi/2  and  y 
in  the  range  of  this  parameter  for  the  solvents  considered  (see  Table  1).  Thus,  one  could  not 
carry  out  a  four  parameter  least  squares  fit  with  In  y^^^,  y,  and  In  Tl  as  independent  variables. 
It  follows  that  the  value  of  g  obtained  in  a  fit  of  eq.  (19)  to  experimental  data  does  not  just 
reflect  the  size-distance  parameter  when  a  is  close  to  unity. 

The  results  of  fitting  eq.  (19)  to  data  for  homogeneous  electron  transfer  are  summarized 

in  Table  2.  Of  the  eight  systems  considered,  the  four  diaminobenzene  systems  can  clearly  be 

considered  to  have  very  small  dependence  on  Tl.  In  fact,  Grampp  and  Jaenicke  [7] 

concluded  that  the  inner  sphere  reorganization  energy  predominated  for  these  systems  and 

analyzed  the  solvent  dependence  on  the  basis  of  eq.  (17).  Although  the  number  of  solvents 

studied  was  not  large,  the  quality  of  the  fit  is  excellent  as  reflected  by  the  high  values  of  the 

regression  coefficient.  In  the  case  of  cobaltacene  (CB)  and  ferrocene,  the  value  of  a  is 

* 

intermediate  lying  close  to  0.5.  Since  these  systems  have  low  values  of  [  10],  this 
result  may  indicate  that  the  reactions  are  weakly  adiabatic.  In  the  case  of  the  cobaltacene 
system,  the  quality  of  the  fit  wai  not  high  when  data  for  all  eleven  solvents  were  considered. 
However,  when  data  for  outlying  points  (BN  and  AC)  were  removed,  the  value  of  the 
regression  coefficient  increased  markedly.  This  could  be  taken  as  evidence  that  the 
coefficient  a  is  changing  with  Tl  since  these  solvents  lie  and  the  ends  of  the  range  of  values 
of  Tl  considered  in  these  experiments  [10].  On  the  basis  of  the  partial  regression 
coefficients,  it  is  clear  that  50  to  60%  of  the  observed  dependence  of  the  logarithm  of  the  rate 
constant  on  solvent  can  be  attributed  to  variation  in  In  Xl-  For  tetracyanoquinodimethane 
(TCQM)  and  decamethylcobaltacene  (DMCB),  the  value  of  a  is  unity  within  experimental 
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error  which  indicates  that  these  systems  can  be  treated  as  adiabatic  with  a  dominating 
contribution  to  the  free  energy  of  activation  from  outer  sphere  reorganization.  In  the  case  of 
the  TCQM  system,  this  conclusion  was  reached  earlier  by  Harrer  et  al.  [23]  who  first 
proposed  the  analysis  of  data  for  such  systems  on  the  basis  of  eq.  (18).  The  results  obtained 
show  clearly  that  the  solvent  contribution  to  the  free  energy  barrier  makes  up  close  to  40%  of 
the  explained  solvent  dependence  so  that  analytical  methods  used  to  separate  the  pre¬ 
exponential  solvent  dependent  factors  from  the  exponential  ones  should  be  designed  to  avoid 
assumptions  about  the  size-distance  parameter. 

Results  of  heterogeneous  kinetic  studies  of  solvent  effects  are  summarized  in  Table  3. 

For  these  systems,  a  wide  variation  in  the  a  parameter  is  also  seen  with  error  estimates  much 

larger  than  for  homogeneous  systems  with  the  same  number  of  solvents.  In  the  case  of  the 

heterogeneous  process,  the  net  change  in  the  rate  constant  with  solvent  is  smaller,  so  that  the 

quality  of  the  fit  is  not  expected  to  be  as  good  as  for  the  corresponding  homogeneous  system 

given  the  same  number  of  solvents  and  the  same  range  of  the  parameters  y.  and  In  Xl-  For  the 

metal  sandwich  compounds,  a  is  unity  within  experimental  error,  and  the  majority  of  the 

solvent  dependence  can  be  attributed  to  the  Xl  dependence  of  the  pre-exponential  factor. 

These  reactions  are  the  fastest  of  those  considered  in  Table  3,  and  were  chosen  to  investigate 

solvent  effects  because  they  possess  low  inner  sphere  reorganization  energies  [16,18].  In 

the  case  of  phenothiazine  and  1,4-diaminobenzene  (DAB)  intermediate  values  of  a  are  found, 

whereas  for  t-nitrobutane  a  has  its  lowest  value.  In  cases  where  data  are  available  for  both 

the  homogeneous  and  heterogeneous  reaction  (CB  and  DAB),  it  is  seen  that  a  increases  in 

going  from  the  homogeneous  reaction  to  the  heterogeneous  one.  This  can  be  attributed  to  a 

*  * 

corresponding  decrease  in  relative  to  AG^^,  or  to  an  increase  in  donor-acceptor 
coupling.  On  the  basis  of  the  present- results,  the  analyses  of  the  heterogeneous  DAB  and 
phenothiazene  data  given  earlier  [15,17,27]  are  incorrect  in  that  it  was  assumed  that  the 
coefficient  a  is  one. 


The  kinetic  data  for  those  systems  in  which  the  parameter  a  can  be  considered  unity 

within  experimental  error  were  reanalyzed  on  the  basis  of  eq.  (18)  in  order  to  determine  the 

size-distance  parameter  Rp  and  the  pre-exponential  parameter,  icKp.  These  results  together 

with  the  molecular  radius  estimated  from  crystallographic  or  density  data  are  summarized  in 

Table  4.  The  kinetic  data  for  DMCB  which  fulfill  the  criteria  for  application  of  eq.  (18)  in 

both  the  homogeneous  and  heterogeneous  reactions  are  also  presented  in  Fig.  1.  It  is 

interesting  that  the  estimates  of  Rp  for  this  system  are  equal  within  experimental  error. 

Thus,  the  slope  of  the  plot  for  the  homogeneous  data  is  approximately  twice  that  for  the 

heterogeneous  data.  The  parameter  Rp  which  is  equal  to  Rna  /  (Rh-a)  in  the  case  of 

homogeneous  reactions,  and  Rea  /  (Re-a)  in  the  case  of  heterogeneous  reactions  is  expected  in 

general  to  be  different  in  the  two  reaction  modes.  The  fact  that  Rp  is  approximately 

independent  of  reaction  mode  in  the  case  of  DMCB  suggests  that  imaging  effects  are 

important  in  the  heterogeneous  reaction.  In  previous  analyses  of  solvent  effects  on 

* 

heterogeneous  electron  transfer  data  [15-18],  has  been  estimated  assuming  that  Rp  =  a, 

that  is  assuming  that  imaging  is  not  important  (Re  =  <»).  On  the  basis  of  the  data 

summarized  in  Table  4,  Rp  is  always  greater  than  the  reactant  radius  a,  so  that  the  assumption 

* 

that  Rp  =  a  can  lead  to  serious  over-estimation  of  AG^^.  Estimation  of  Rp  is  not  simple  for 

* 

polyatomic  molecules  under  any  circumstances.  This  aspect  of  the  variation  in  AG^^  with 

solvent  has  been  discussed  by  Jaenicke  and  coworkers  [7-9]  who  have  demonstrated  the 

importance  of  considering  carefully  molecular  shape  in  estimating  the  size-distance 

parameter.  This  feature  of  electron  transfer  theory  becomes  more  important  with  increase  in 

the  number  of  atoms  in  the  reactant  since  the  mutual  orientation  of  the  reactants  in  the 

transition  state  must  be  considered,  electron  transfer  often  being  associated  with  a  specific 

location  in  the  molecule.  Nevertheless,  in  spite  of  the  fact  that  the  tl  solvent  dependence 

plays  a  dominant  role  in  determining  the  net  solvent  effect,  the  present  results  provide  very 

* 

strong  confirmation  of  the  Marcus  model  for  AG,... 
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It  is  also  interesting  to  compare  the  values  of  the  pre-exponential  parameter,  icKp,  with 

that  predicted  on  the  basis  of  the  encounter  pre-equilibrium  model  [25,30].  Assuming  that 

Rh  =  2a  and  5r  =  60  pm  [30],  icKp  for  homogeneous  adiabatic  reactions  varies  between  0.16 

and  0.45  dm^  mol'^  for  variation  in  the  reactant  radius  between  0.3  and  0.5  nm  (see  eq.  (6)). 

The  value  of  icKp  obtained  for  TCQM,  namely,  0.6  dm^  mol'^  is  slightly  higher  than 

expected  on  the  basis  of  the  parameters  listed  in  Table  4,  but  is  of  the  correct  order  of 

magnitude.  In  the  case  of  DMCB,  the  estimate  of  icKp  is  higher  than  expected  by  a  factor  of 

30.  This  can  be  partially  attributed  to  the  poor  quality  of  the  fit  which  may  lead  to  a  low 

value  of  the  size-distance  parameter  and  thus  a  high  value  of  icKp.  If  one  forces  the  size- 

distance  parameter  to  be  0.96  nm  (Rp  =  2a),  the  estimate  of  icKp  drops  to  1.3  dm^  mol'^  a 

result  which  is  closer  to  the  range  of  values  predicted  by  theory.  In  fact,  Nielson  et  al  [10] 

achieved  agreement  between  experimental  and  theoretically  predicted  values  of  this  parameter 

when  it  was  assumed  that  Rp  =  1 ,3  nm.  In  order  to  improve  the  quality  of  the  fit,  data  are 

needed  in  more  soivents  with  an  improvement  in  the  precision  of  the  estimates  of  Tl.  The 

latter  point  is  discussed  in  more  detail  below.  In  the  case  of  the  heterogeneous  data,  the 

experimental  values  of  xKp  all  fall  below  the  predicted  value  of  60  pm  for  adiabatic  reactions 

[30].  In  analyzing  these  data,  it  has  been  customary  in  previous  work  [15-18]  to  assume 

that  the  size  distance  parameter  Rp  is  equal  to  the  molecular  radius  a,  and  to  assess  the  solvent 

* 

dependence  using  values  of  estimated  accordingly  (see  eq.  3  ).  In  the  case  of  the 
cobaltacene  system  for  which  the  most  data  are  available,  the  estimate  of  icKp  assuming  Rp  = 
0.37  nm  is  380  pm,  a  result  which  is  an  order  of  magnitude  higher  than  that  predicted  by 
theory.  A  fact  that  should  be  considered  in  assessing  the  heterogeneous  results  is  that 
double  layer  effects  have  been  ignored  in  analyzing  the  data  for  most  systems  [15-18]. 

Thus,  the  true  or  double  layer  corrected  pre-exponential  factor  may  be  higher  or  lower  than 
that  given  in  Table  4  depending  on  whether  the  electrode's  field  accelerates  or  decelerates  the 
reaction. 
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Kinetic  data  for  which  a  can  be  assumed  zero  were  reanalyzed  according  to  eq.  (17),  the 
results  being  summarized  in  Table  5.  The  four  systems  which  fulfill  the  criteria  for  application 
of  eq.  (17)  are  the  homogeneous  electron  transfer  reactions  involving  diaminobenzenes  studied 
by  Grampp  and  Jaenicke  [7].  As  reported  earlier  [7],  the  size-distance  parameter  for  these 
systems  is  much  larger  than  twice  the  molecular  radius,  that  is,  the  value  expected  for  spherical 
reactants  at  contact.  This  feature  can  be  attributed  to  the  non-spherical  nature  of  the  reactant 
and  has  been  discussed  in  detail  earlier  [7].  What  is  more  interesting  here  is  the  relative  values 

Ik  Ik 

of  AG  and  AG  and  the  reason  why  these  systems  exhibit  little  or  no  dependence  of  the 
pre-exponential  factor  on  tL-  On  the  basis  of  the  size-distance  parameters  recorded  and 

Ik 

choosing  an  average  value  of  the  permittivity  parameter  y  (0.4),  the  average  value  of  AG  is 

* 

approximately  twice  that  of  AG Thus,  the  Ovchinnikova  criterion  [3,13]  is  not  met  for 

most  of  the  solvents  involved.  In  going  from  the  homogeneous  to  the  heterogeneous  system, 

one  finds  an  increase  in  the  Tl  contribution  to  the  solvent  effect  in  the  case  of  DAB.  This 

*  ★ 

cannot  be  attributed  to  an  increase  in  AG  with  respect  to  AG  unless  there  is  a 
corresponding  decrease  in  the  size-distance  parameter  by  a  factor  of  two,  or  unless 

Ik 

orientational  requirements  manifested  in  the  respective  values  of  AG  only  are  very  different 
for  the  electrode  reaction  with  respect  to  the  homogeneous  one. 

In  assessing  the  above  results  one  should  keep  in  mind  that  the  parameter  xl.  which 
accounts  for  the  major  portion  of  the  solvent  effect,  is  often  imprecisely  known  and 
depends  on  the  nature  of  the  indifferent  electrolyte  [35].  Although  the  Debye  relaxation 
time  to  can  be  determined  with  reasonable  precision,  the  infinite  frequency  dielectric 
constant  £«  has  often  been  estimated  on  the  basis  of  dielectric  relaxation  data  obtained  in 
too  low  a  frequency  range,  long  extrapolations  to  frequencies  up  to  ~1  THz  being  required. 

As  an  example,  we  cite  the  data  for  propylene  carbonate  [41-44],  an  aptotic  solvent  with  a 
high  dielectric  constant  and  long  relaxation  time.  The  available  data  which  are  summarized 
in  Table  6  demonstrate  clearly  that  the  estimates  of  Xl  vary  over  an  unacceptably  wide  range 
and  that  the  major  source  of  imprecision  is  the  high  frequency  dielectric  constant  £oo.  As 
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the  Debye  relaxation  time  decreases,  the  problem  becomes  worse  because  measurements  at 
higher  frequencies  are  required  to  determine  Eoo.  In  fact,  for  solvents  with  fast  relaxation 
times  such  as  acetonitrile,  the  values  of  £«  are  only  estimates,  no  experimental  values  being 
available.  The  problem  addressed  here  is  apparent  in  the  discussion  of  other  experimental 
results,  especially  those  related  to  time  resolution  studies  of  fluorescence  phenomena  in 
different  solvents  [45,46].  It  is  clear  that  more  data  for  solvent  relaxation  processes  are 
required  in  a  higher  frequency  range  both  in  the  presence  and  absence  of  electrolytes  before 
the  questions  raised  here  can  be  resolved.  In  the  analyses  presented  here,  when  more  than 
one  estimate  of  Tl  is  available  from  the  literature  the  value  chosen  was  that  which  resulted  in 
the  best  fit  of  eq.  (16)  to  the  kinetic  data  (see  Table  1). 

On  examining  the  literature  dealing  with  experimental  data  for  solvent  effects  on 
electron  transfer  reactions  [6-19],  it  is  readily  apparent  that  conclusions  qu’te  different  than 
those  reached  here  have  been  reported  in  many  cases.  This  situation  has  arisen  because  it 
has  been  common  practice  to  assume  a  certain  solvent  dependence  for  the  free  energy 
barrier  in  the  rate  constant  and  then  examine  the  solvent  dependence  of  the  pre-exponential 
factor.  We  emphasize  that  the  procedure  described  here  allows  one  to  separate  the 
dependences  on  the  permittivity  parameter  y  and  the  longitudinal  relaxation  time  Tl  subject 
to  the  number  and  quality  of  the  corresponding  data,  and  thus  reach  statistically  more  sound 
conclusions  regarding  roles  of  these  two  parameters  in  determining  the  observed  rate 
constant.  For  instance,  Nielson  et  al.  [10]  concluded  that  kinetic  data  for  DMCB  in  the 
homogeneous  mode  have  a  fractional  dependence  on  tl  (0.7  >  a  >  0.85)  whereas  the 
present  analysis  shows  that  a  is  much  closer  to  unity.  This  disagreement  is  almost 
cenainly  due  to  the  fact  that  these  authors  assumed  the  dependence  of  the  kinetic  data  on  y 
to  obtain  the  dependence  on  xl.  Since  the  dependence  on  y  makes  up  close  to  40%  of  the 
total  explained  variation,  their  analysis  is  clearly  subject  to  doubt. 

In  conclusion,  the  value  of  the  present  analysis  is  clearly  dependent  on  the  validity  of 
eq.  (12).  We  have  already  pointed  out  above  that  in  some  cases  the  parameter  a  may 


depend  on  solvent  nature  for  the  range  of  solvents  considered  in  a  given  analysis.  In  this 
case,  our  analysis  gives  a  good  indication  of  the  relative  importance  of  the  two  parameters 
in  determining  the  solvent  effect  but  the  derived  value  of  a  must  be  considered  an  average 
value.  There  are  other  problems  of  a  more  subtle  namre  that  could  be  considered,  one  of 
them  being  the  fact  that  Re  is  expected  to  change  with  solvent  nature  for  heterogeneous 
reactions  due  to  a  corresponding  variation  in  the  distance  of  closest  approach  of  the  reactant 
to  the  electrode  with  solvent  size  [11].  However,  the  most  important  problem  is  the  poor 
quality  of  the  tl  data  for  many  of  the  solvents  considered  in  current  experimental  studies. 
Further  experimental  work  investigating  solvent  relaxation  behaviour  at  high  frequencies 
would  do  much  to  improve  our  understanding  of  the  solvent's  role  in  the  kinetics  of 
reactions  in  solution. 
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Table  1.  Solvent  Parameters  Relevant  to  Estimation  of  the  Electron  Transfer  Rate 


Constant  (Eq.  19) 


Solvent 

Longitudinal 

Permittivity 

Relaxation  Time^ 

Parameter^ 

ps 

Y 

Acetone  (AC) 

0.3 

0.495 

Acetonitrile  (AN) 

0.2 

0.529 

Benzonitrile  (BN) 

5.8 

0.390 

Chloroform  (CF) 

2.4 

0.276 

Dichloroethane  (DCE) 

1.6 

0.384 

Dichloromethane  (DCM) 

0.9 

0.382 

Dimethoxyethane  (DMXE) 

0.8 

0.371 

Dimethylacetamide  (DMA) 

1.5 

0.459 

Dimethylformamide  (DMF) 

1.1 

0.463 

Dimethylsulfoxide  (DMSO) 

2.1 

0.437 

Hexamethylphosphoramide  (HMPA)  8.8 

0.438 

Nitrobenzene  (NB) 

5.3 

0.387 

Nitromethane  (NM) 

0.2 

0.498 

Propylene  Carbonate  (PC) 

1.7 

0.480 

Pyridine  (PY) 

1.3 

0.359 

Tetrahydrofiiran  (THF) 

1.7 

0.388 

Tetramethylurea  (TMU) 

6.0 

0.433 

3  Defined  in  eq.  (10). 
Defined  in  eq.  (14). 


Table  4.  Kinetic  Parameters  for  Adiabatic  Electron  Transfer  Reactions  for  which  AG;  «  AG  (Equation  1 8). 
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a  Calculated  from  density  or  crystallographic  data  assuming  that  the  volume  per  molecules  is  47ca^/3. 
The  units  of  Kp  are  dm^  mol‘*  for  homogeneous  reactions  and  pm  for  heterogeneous  reactions. 


Table  5.  Kinetic  Parameters  for  Electron  Transfer  Reactions  Analyzed  According  to  Equation  17. 
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**  Based  on  estimates  made  by  Grampp  and  Jaenicke  [7]  and  revised  by  Kapturkiewicz  and  Jaenicke  [40]. 


Table  6.  Dielectric  Relaxation  Parameters  for  Propylene  Carbonate  at  25 
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Legends  for  Figures 
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Figure  1 .  Plot  of  the  solvent  corrected  kinetic  parameter  In  (k  against  the 

permittivity  parameter  y  using  kinetic  data  for  homogeneous  electron  transfer 
(O)^®  and  heterogeneous  electron  transfer  (A)^^  for  the  decamethylcobaltacene  - 
decamethylcobaltacenium  system.  The  left  hand  ordinate  scale  applies  to  the 
homogeneous  data  and  the  right  hand  scale,  to  the  heterogeneous  data. 
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